ASYMPTOTIC EXPANSIONS AND EXTREMALS FOR THE CRITICAL 
SOBOLEV AND GAGLIARDO-NIRENBERG INEQUALITIES ON A TORUS 



MICHELE BARTUCCELLI, JONATHAN DEANE AND SERGEY ZELIK 

Abstract. We give a comprehensive study of interpolation inequalities for periodic functions 
with zero mean, including the existence of and the asymptotic expansions for the extremals, 
best constants, various remainder terms, etc. Most attention is paid to the critical (logarithmic) 
Sobolev inequality in the two-dimensional case, although a number of results concerning the best 
^ ' constants in the algebraic case and different space dimensions are also obtained. 
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1. Introduction 
We study the following critical Sobolev inequality (see [5l [15]): 

(1-1) Ml^iu) < Mmhu) (ci log ^J^^ + c)i 

in the particular case when is a two dimensional torus. This inequality, which can be formally 
considered as a limit case (/ — )■ 1, n = 2, d=2) of the algebraic inequality of the Gagliardo- 
Nirenberg type: 

(1.2) Wuh^^n) < Cf,(/,n)||(-A)-'/\||^.(^)||(-A)-"/\||i^(%, 

^ ^ (n - d/2) ^ ^>^/2>/, nccR'^ 
n — I 

is known to be very useful in many problems related to partial differential equations and mathe- 
matical physics. For instance, it is used for obtaining best known upper bounds for the attractor 
dimension of the Navier-Stokes system on a 2D torus (see e.g. |29J), for proving the uniqueness 
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of weak solutions for von Karman-type equations arising in elasticity (see and references 
therein) as well as for the so-called hyperbolic relaxation of the 2D Cahn-Hilliard equation (see 
[14| ) or the 2D Klein-Gordon equation with exponential nonlinearity (see [I?]). We mention 
also that a slightly different logarithmic inequality is used at a crucial point in the proof of the 
global existence of strong solutions of the 2D Euler equations (see [33j). 

Note that nowadays most classical inequalities of Gagliardo-Nirenberg type can be relatively 
easily verified using interpolation theory (see e.g., [30]). However, the best constants in those 
inequalities as well as the existence and the analytic structure of the extremals is a much more 
delicate and interesting question, which is far from being completely understood despite per- 
sistent interest in the problem and the many interesting results obtained during the last 50 
years; see [D 13 El El [IHl [III IIS III |23l [28l EZl EI] and references therein. The most studied 
is, of course, the case of the whole space Q = W^, more or less complete results are available in 
two cases: where the inequality does not contain derivatives of order higher than one or in the 
Hilbert case. In the first case the rearrangement technique works and reduces the problem to the 
one-dimensional case and in the second case one can use the Parsevsal equality. In particular, 
as proved in [19], the best constant in (11. 2p for the case = R*^ is 

(1.3) cjjd(/,n) = 



1 X 1/2 

7ruj{d) f 1 \ 2("-i) \ 



^(2vr)'isin2^ \{d - 2lY-^^i2n - d)^^-^ 

where u}{d) = p^'^^^ is the surface area of the (d — l)-dimensional sphere. In addition, the 

extremal function G (-A)^"/^L^(M'^) fl (-A)^'/^L^(IR'^) exists and is unique up to a shift 
and scaling n*(x) — )• an*(/3x — xq), a,/3 G M, xq G W^; n* is given by 

(1-4) Mx) = j^lj^^;^e-^dx. 

The situation becomes more complicated in the case where is a bounded domain of W^, even 
for the algebraic inequality (|1.2p with Hilbert norms on the right-hand side. To the best of our 
knowledge, two different scenarios are possible here. In the first case, the sharp constant cn{l, n) 
coincides with Cj^d{l,n), but in contrast to the case of W^, there are no exact extremals and the 
approximative extremals can be constructed by the proper scaling and cutting of the function 
(jl.4p . This case is realized, for instance, if Dirichlet boundary conditions are posed; or if $7 = 
is a circle (periodic boundary conditions) and n = 0; or if = S'^ is a higher dimensional sphere 
(d = 2 and A is a Laplace-Beltrami operator) with n = and / < 7. See jl8t [T9] for details. 
In the present paper, we show that it is also true for the tori = and = if Z = and 
n is not too large — see Section El In addition, in that case, inequality ()1.2p can be improved 
by adding an extra lower order term in the spirit of Brezis and Lieb (see 0). In particular, as 
shown in Section El the following inequality 

(1-5) ll^llioo(']r2) < -||'u||2,2(x2) • ||Aa;n||i2(T2) - ^;^ll«lli2('ir2) 

holds for all 27r x 27r-periodic functions with zero mean. However, even for this improved 
inequality the exact extremal functions do not exist, and further improvements can be obtained. 

In the second case, the sharp constant in (11.21) is strictly larger than the analogous constant 
in M"^: 

(1.6) cn{l,n) > c^d{l,n) 

and there is/are exact extremal function(s) for ()1.2p in i7"(il). In particular, this holds for 
= §2 with / = 0, n > 8, see [19] (see also [18] for the analogous effect for the slightly different 
inequality in the one-dimensional case). In that case, the constant CQ(/,n) can be found only 
numerically as a root of a transcendental equation. 
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It was conjectured by Ilyin that the analogous effect holds on multi-dimensional tori = T"^, 
d > 1 and ^ = 0. In the present paper, we verify that this conjecture is indeed true and that 
the inequality (|1.6p holds for 17 = (with zero mean) for / = and n = 10. In addition, we 
establish the following 3D analogue of (jl.Sp : 



II ||2 ^ \/2\/3 ,, 1/2 11 ,,3/2 ^11 ||2 

) IPIIl°°{t3) - ^ \m\L'2(T3')\\^'^\\L^{Ti) ~ ^ II'"IIl2{t3)' 

where the sharp constant still coincides with the analogous constant in the whole space 

R^, but nevertheless (II. 7p possesses an exact extremal function and the best value for the second 
constant K can be found only numerically {K ~ ^^J^ , see Section [5]). 

Let us now return to the limit logarithmic inequality (jl.ip . This case looks more difficult than 
the algebraic one, in particular, since it is a priori not clear whether or not the transcendental 

||„||2 

function 5 — )■ Ci log (5 + C2 (5 := ,, nf ^"^ ) on the right-hand side of (jl.ip is optimal. Indeed, a 
detailed study of the slightly different logarithmic inequality 

V nn^Q) J 

where the H^-norm. is replaced by the Holder norm with a G (0, 1) is given in recent papers 
[21 [16] for the case where Q is a unit ball and the function u satisfies the Dirichlet boundary 
conditions (see also [HI [25l [26]). As shown there, Ci > and, in order to be able to take 
C'l = an extra double logarithmic corrector [5 log log (5) is required. Thus, based on that 
result and on the interpolation inequality 

one may expect the following improved version of (jl.ip 



(1.9) 11^112^(1^2) < ^||V«||^2(Tr2) 



jroved version 01 (|l.ip 

,(log^*ll!i + loJl + l„g^|!<I!i]+L],L>0 

\ II vu||^2(Tr2) y II^"IIl2(q) J J 

T X 27r-periodic functions u with zero mean. Note that tl 
ied on the reducine the problem to the radially symmetric cai 



to be optimal for the case of 27r x 27r-periodic functions u with zero mean. Note that the 
analysis presented in ^ [16] is based on the reducing the problem to the radially symmetric case 
via the rearrangement technique, and use of the Dirichlet boundary conditions is essential, so it 
is not clear how to extend it — either to the case of the torus or to the case of the H^-norm. 
Nevertheless, as we will see below, inequality ()1.9p is true for the properly chosen constant L 
(which can be found numerically as a solution of a transcendental equation: L ~ 2.15627). In 
addition, there exist exact extremal functions for this inequality; see Section [3l 

The main aim of the present paper is to introduce a general scheme which allows the analysis of 
inequalities (|l.ip . (|1.2p and (jl.Op at least on tori, from a unified point of view, and to illustrate it 
in the most complicated logarithmic case (although nontrivial applications to the algebraic case 
will be also considered). One of the important features of our approach is that, in contrast to, 
say, |16[ [To] (and similarly to [3]), the concrete form of the right-hand sides in those inequalities 
is not a priori postulated, but appears a posteriori as a result of computations. Indeed, instead 
of (II. 9p . we consider the following variational problem with constraints: 

llniP r IIA-ulP, 

(1.10) , ' ^ max, u G ^^(T^), / u{x) dx = 0, I ^ = 5 

||Vn|||^2 Jt2 11^^112^2 
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and prove that, for every 6 > 0, this problem has a unique (up to shifts, scahng and alternation 
of sign) solution u^{x), 

(1.11) u^{x)= k^i + fik^y ^'^=^? + ^2 

(compare with (jl.4l) ) and the parameter fi can be found, in a unique way, as a solution of the 
equation 

Let us denote the maximum in (jl.lOp by @{S); as we will see, is a real analytic function of 6. 
Then, the following inequality holds: 

/\\Au |p2 

(1-13) ll^^lli-(T2) < l|Vn|||2(Tr2)e f 

and by definition is the least possible function in this inequality. Thus, inequality (jl.lSp 
can be considered as an optimal version of (jl.ip and (jl.9p . However, inequality (|1.13p is not 
convenient for applications since the function Q is given in a very implicit form through the 
lattice sums (II. lip which, to the best of our knowledge, cannot be expressed in closed form 
through the elementary functions (in contrast to the case of inequality (jl.Sp in a unit ball, see 
[3j) and, in addition, direct numerical computation of them is not easy especially for large 5 
(small n) due to very slow rate of convergence. 

In order to overcome this problem, we have found the asymptotic expansions for the function 
G((5) as (5 — )• oo. Namely, we have proved that the function Q(5) coincides up to exponentially 
small terms (of order 0{e~'^'^^^ ^)) with the function 0o(5) given by the following parametric 
expression: 

^ ^ ^TT^ vrlogi +;S-7r + 2;u' vrlogi +/3-7r + 2//' 

where /3 := vr (27 + 21og2 + Slogvr — 41ogr(l/4)), 7 is the Euler constant and T(z) is the Euler 
gamma function. In particular, 

Go(5) = ^ log 5 + ^ log log 6 + + 05^00(1), 

which justifies inequality (jl.9p and shows that the constant L > L^o := In practice, the 

numerics shows that L > Lopt > L^o — see Section [3l 

In addition, combining the analytic asymptotic expansions for Q{8) with numerical simulation 
for relatively small 5, we show that 

(1.15) e{5) < Qo{5) 

for all 6 >!. Thus, the much simpler function Gq can be used instead of G in the right-hand side 
of (jl.lSp . Actually, Go gives a reasonable approximation to G for all values of 6. For instance, 
for <5 = 1,2 and 4, respectively, we have G (resp. Go) = 0.10134(0.17797), 0.26651(0.26660) and 
0.35112(0.35112). 

The paper is organised as follows. The proof of the existence of the conditional extremals for 
problem (jl.lOp as well as analytical formulae for them in terms of the lattice sums, are given in 
Section [2j 

The key asymptotic expansions for the lattice sums involving the parametric expression for 
G((5), as well as for the extremals u^{x), are presented in Section[3l Based on these expansions, 
we check the validity of inequality ()1.9p as well as the estimate ()1.15p . 
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The elementary approaches to the logarithmic inequality (jl.ip are analyzed in Section [H 
Actually, there are at least two known ways to prove this inequality without studying the 
corresponding extremal problem: one of them is based on the embedding H^'^'^ C C with 
further optimization of the exponent e > (see e.g. [2]), and the other more classical one 
(which has been factually used in the original paper [5]) splits the function u into lower and 
higher Fourier modes and estimates them via the and H"^ -norms respectively. Based on the 
above asymptotic analysis, we show that the second method is preferable and allows to find the 
correct expressions for the two leading terms in the asymptotic expansions of the function Q. 

The application of our approach to the simpler algebraic case ()1.2p with / = and arbitrary 
space dimension d is considered in Section [3 We establish here the following improved version 
of (O) 



(1-16) < c,(n)h||5^;''/"||(-A,r/2n||g" - K,{n)\\u\\l2, 

where Cd(n) = c^d{0,n) for all n G N such that 2n — d > and the constant Kd{n) may be 
either positive or negative. We prove that in the one dimensional case this constant is strictly 
positive, but it may be either positive or negative in the multi-dimensional case, in depending on 
n. We present also combined analytical/numerical results for the constants Kd{n) for d and n 
not large. In particular, inequalities (jl.Sp . (|1.7p mentioned above, as well as the ID inequalities 

2 II /II -''ii ii2 ii2 1 1 3/2 II //II 1/2 2 2 

||u|b(Ti) ^ ||l2 - -\\u\\l2, \\u\\c(T^) ^ :^II^IIl2 \\u 11^2 - ^\Ml^, 

are verified there. 

The large n limit of the inequality (jl.lGp is studied in Section [6l The results of this section 
clarify the nature of oscillations in the analog of the function (5 — )• 0((5) for that inequality and 
show the principal difference between the ID case where the regular oscillations occur (after the 
proper scaling) and the multi-dimensional case where the oscillations are irregular due to some 
number theoretic reasons. 

Finally, the computation of the integration constant /3 is given in the Appendix. 

Acknowledgement. The authors would like to thank A. Ilyin for many stimulating discus- 
sions and comments. 

2. Conditional extremals: existence, uniqueness and analytical expressions 

This section is devoted to the study of the maximisation problem (jl.lOp which we rewrite in 
the following equivalent form: 

(2.1) ||n||^(Ti.2) ^ sup, ueH'^{T'^),j u{x)dx = 0, ||An|||2 = <5, ||Vn|||2 = 1. 

In addition, we note that problem (j2.ip is invariant with respect to translations u{x) — )• u{x + h) 
and alternation u{x) — )■ —u{x). Thus, without loss of generality, we may assume that ||ii||c(T2) = 
u(0) > and so reduce problem ()2.ip to the following one: 

(2.2) tt(0)^sup, ueH^{T^), / u{x)dx = 0, ||An||^2 = -5, ||Vn||^2 = 1. 
Thus, the function in (11.13P can be defined as follows: 

(2.3) 9(5) :=sup|n(0)2, ueH^{T^),j u{x)dx = 0, ||An|||2 = -5, ||Vu||^2 = l|. 



It is however more convenient to rewrite problem (12. 2p and (|2.3p in Fourier space by expanding 



where Yl' means the sum over the lattice k ^1? except A; = 0. Using the Parseval equality, we 
transform (j2.2p to 

(2.5) ^^'ufc^sup, , {k^ f\uk\^ = 5, k^\uk\^ = I. 

Finally, we observe that, without loss of generality, we may assume that all Uk in (I2.5p are real 
and nonnegative. 



Lemma 2.1. For every S >1 there exists an extremal function (maximiser) for problem (12.2p 
(or equivalently, for problem (12.5P ). 

Proof. Let Un{x), Un{0) > be a maximising sequence for problem ()2.2p such that 

e{6) = lim Un{Of. 

Such a sequence exists if and only if 5 > 1, since under that condition the set of functions 
u £ ff^(T^) for which the constraints of ()2.2p are satisfied is not empty. Clearly, n„ is bounded 
in and consequently, without loss of generality, we may assume that n„ — )• u* weakly in 
(and strongly in C(T^) and in H^). We claim that u* is the desired maximiser. Clearly, 

(2.6) 9(5) = u*{Of, \\Vu*\\l2 = 1, \\Au4l2 < 6. 

Thus, we only need to check that the last inequality is in a fact equality. Assume that it is not 
true and ||A2,.u*||^2 = Sq < 6. Let us fix ko G such that Uk^ > take any small e > and 
A^ > 1^0 1 and consider the perturbed function 



where /3 = l3{e,N) > is chosen in such way that ||Vne^Ar||^2 = 1. Using the fact that 
Y' \k\'2 iog\\k\+i) ^ °° ^^^^ > 0' '^^^ easily show that there are positive constants eo 
and / independent of A^ such that 

(2.7) l3{e,N) < le, Ve < eo 

and all A^. On the other hand, using (j2.7p and the fact that iog(|fc|+i) = 00, we see that 

there exists Nq independent of e such that 

(2.8) Ue,N{0) > u^{0), yN>No,e<eo. 
Finally, since 

lini ||Aa;M£ Ar||^2 = Sq < 5, lim || A^.^^ Ar||i2 = oo, 

we may find A'^* > A''o and e^, < eo such that ||Aa;U£,^Ar, ||^2 = <5- This, together with ()2.8p . shows 
that 

e{5) > n,(0)2 

which contradicts our choice of function u=k (x) (see ()2.6p ) and finishes the proof of the lemma. □ 

Remark 2.2. In particular, the above arguments show that the function 6 — )• Q{5) is strictly 
increasing. 

We are now ready to state the main result of this section, which gives the existence and unique- 
ness for the extreme functions of ()2.ip . These functions will be further referred as conditional 
extremals for the logarithmic Sobolev inequality considered. 
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Theorem 2.3. For every fixed 5 > 0, variational problem (jl.lOp has a unique (up to translations, 
scalings and alternation) solution 

where fi = fi{5) £ (— oo, —1] U (0, oo] is determined as the unique solution of the equation 

(2.10) Fi,) := 5^ = 

Thus, the desired function Q{5) possesses the following parametric representation: 

^ 1 |fcP(l+HfcP)) . ^' (l+Mfc^)^ 

(2.11) QM-=X^- v^/ 1 — ' ^i^") ■= ^^ r-^i 

and fi G (— oo, —1] U (0, oo]. 

Proof. Instead of (jl.lOp . we will consider the equivalent problem (j2.5p . The extremals of that 
problem can be easily found using Lagrange multipliers. Introducing the Lagrange function 

C{u) := ^Y^uk + Ai \k\\l + A2 \k\\l, Ai,A2e M, 

differentiating it with respect to and using the necessary condition ■^C{u) = for extremals, 
we find the following extremals: 

^^•^^^ " ""^^'^i'^^ = 47r|A;|2(^i+yl2|fc|2)' 

where, as usual, the multipliers Ai and A2 should be chosen to satisfy the constraints. Since we 
already know (from Lemma l2.ip that the maximiser us{x) exists, its Fourier coefficients should 
satisfy p.l2p for some Ai and A2- Moreover, taking into the account the fact that the initial 
variational problem is scaling invariant, we may get rid of one of the multipliers Ai and A2 by 
introducing /x = A2/A1. We will then end up with the one-parameter family of extremals (12. 9p 
depending on /i (the case j4i = is not lost and will correspond below to /i = 00). Of course, 
the parameter /i should be chosen to satisfy the constraints, namely, || An^||22/||Vu^||^2 = 5- 
This gives equation (|2.10p . and representation (j2.1ip follows immediately from the definition of 

Thus, we only need to verify that the solution of F{ij) = 5\s unique. To this end, we first recall 
that all the Fourier coefficients of the conditional maximiser (s) should be either non- negative 
or non-positive. This, together with the formula (j2.9p . gives the condition that /x G (—00, —1] 
which corresponds to all negative coefficients and ^ G (0, 00] which corresponds to all positive 
ones. Thus, only the values S (—00, —1] U (0, 00) may correspond to the true maximisers and 
we need not consider the case /x S (—1,0). The following Lemma gives the uniqueness of a 
solution of (I2.10p in this domain of /i. 

Lemma 2.4. The function F{e) := F{e^^) is continuous (in fact real analytic), strictly increas- 
ing on [— 1,00) and satisfies 

(2.13) F(-l) = 1, lim F{e) = +00 

e— >+oo 

Therefore, the solution of F{fi) = 6, fi £ (—00, —1] U (0, 00), exists and is unique for all 6 >1. 
Proof. The function F{e) is 

Fie) 



S |A;|2(l+e|fc|2) 



and, differentiating this with respect to e, we have 

,1 ( 1 , 1 ^ v^/ v^/ 1 fc^-;2 



F'{e) 



The double-double sum in the numerator can be rearranged to contain only positive terms. 
Indeed, putting together the terms corresponding to the indices {l,k) and {k,l), we see that 



Thus, we only need to verify (j2.13p . The first assertion is obvious since both nominator and 
denominator in the definition of F have simple poles at e = — 1 with the same residue. The 
second limit is a bit more difficult, but we do not want to prove it here since the detailed 
analysis of the asymptotic behavior of -F as — >• will be given in the next section. Lemma [2.41 
is proved. □ 

Thus, due to the proven uniqueness of a solution of ()2.10p . the conditional maximizer u^{x) = 
for the variational problem (jl.lOp is also unique and Theorem 12.31 is proved. □ 

Remark 2.5. It is not difficult to see that the extremals u^{x) defined by (j2.9p satisfy the 
following boundary value problem: 

(2.14) A^(l - ^A^)n^ = -4^2^(x) + 1 

endowed by the periodic boundary conditions (here 5{x) is a standard Dirac delta- function) 
and, therefore, are closely related to fundamental solutions for this family of 4th order elliptic 
differential operators. It can also be derived by applying the method of Lagrange multipliers 
directly to problem ()2.2p (without passing to Fourier space). We will use this fact below in order 
to find good asymptotic expansions for u^{x). 

3. Asymptotic expansions 

In this section, we deduce the asymptotic expansions up to exponential order for the lattice 
sums used in Theorem 12. 3| which are crucial for our approach. Namely, we will analyse the 
asymptotic behaviour of the following three sums: 

/(/^) = E'^2(r^' 5(/^) = E' ^2(1+ ^^2)2' MM) = E'(rTW 

as ^ — >• 0. Actually, these sums are closely related to each other: 

(3.2) /'(^) = -/i(^), g{^,) = /(/.) - /./i(^) 

and therefore, up to a nontrivial integration constant, we only need to study the simplest func- 
tion, /l(/U). 

Lemma 3.1. The function h{fi) possesses the following asymptotic expansion: 

(3.3) hi^i) = - - 1 + 4^V"'/^e-2-/v^ (1 + 0^(1)) 
as ^ —)• 00 . 

Proof. The derivation of the expansion is based on the Poisson summation formula. Using the 
fact that 

1 



(3.4) Yl Mk) = E M^'^k), ^^{2 
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(Pfj,{z) = Lpi{p}l'^z) and ^(^) = n ^^{fJ. ^^"^0 together with the fact that ipi is analytic in a 
strip |Im2:| < 1, we conclude that ^(^) is exponentially decaying: 



(61 <ae-(^-")l"l,Ve>0, 



Thus, if we need the asymptotic expansion of the left-hand side of (|3.4p up to exponential order, 
only the term with A: = is needed in the right-hand side and, therefore, replacing the sum by 
the corresponding integral gives an exponentially sharp approximation to the lattice sum: 



(3.5) Mm) = E' (TtW ^ ^'^''^ " ^ 



iR2 (1 -h/i|2;|^)^ /i 



for arbitrary e > 0. However, if we need the leading exponentially small term in expansions 
like (|3.3p . we need to look at four more terms on the right-hand side of (j3.4p . namely, the 
terms corresponding to A; = (0, 1), (0, — 1), (1, 0), (— 1, 0) (other terms will decay faster than 
e^~'^~^^^^'^/ ) . To this end, we need to compute the 2D Fourier transform of the radially 
symmetric function 9J^(|2;|). This can be done, for instance, by noting that that Fourier transform 
is a radially symmetric fundamental solution of the squared Helmholtz operator 

(1 - fiA^fu = 6{x). 

The radially symmetric solution of this equation can be explicitly written in terms of Bessel 
functions: 

(3.6) Rm) ■■= m) = 7T\c\^-'^'K,m/vT^) 

where Ki is standard Bessel X-function of order 1 (see, e.g., [32j). Thus, we only need to find 
the leading term in R{2'k^~^/'^) as /i — oo, which can easily be done by using known expansions 
for the Bessel functions (see [32j): 

R{2^^x''") = vrV-'/'e-^'^/^ (1 + o^(l)) . 

Taking into account the fact that there are four identical terms on the right-hand side of ()3.4p 
which correspond to \k\ = 1, we arrive at (j3.3p and finish the proof of the lemma. □ 

As a next step, we need to derive analogous expansions for f{fi) and g{fj,). However, the trick 
with the Poisson summation formula is not directly applicable here since the corresponding 
function ip will have singularity at x = and, as we will see below, this leads to an extra 
residual- type term in the expansions. Instead, we will use relations ()3.2p in order to find the 
expansions for / and g up to an integration constant. 

Corollary 3.2. The functions f{fJ-) and g{p) possess the following asymptotic expansions as 

(3.7) /(/i) = 7rlog- + ^ + /3 - 47r//i/4e-2-/v^(l + o^(l)) 
and 

(3.8) g{fi) = vr log - + 2^ + /3 - vr - 47rV~'/'e-2-/v^(l + o^(l)), 

A* 

where the integration constant f3 = 7r{2j + 2 log 2 -)- 31og7r — 4 log r( 1/4)). 
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Proof. Up to the integration constant /3, expansions (|3.7p and (j3.8p are straightforward corol- 
laries of p.Sp and (|3.2p . so we only mention here the explicit expression for the integral of the 
leading exponential term in (jS.Sp with respect to fi: 

POO 

/ 47r2a;-5/4e-2-/v^ = 4^2^erf (^/2^/i-i/4), 
J fi 

where erf(a;) is the usual probability integral. Then, using the well-known expansions for the 
erf(z) near z = oo, we find the leading exponential term in (j3.7p (and ()3.8p follows immediately 
from the second formula of ()3.2p ). 

Thus, we only need to find the integration constant /3. This, however, is a much more delicate 
problem and the arguments above do not indicate how to compute it. The derivation, based on 
the Hardy formula for the 2D analogue of the Riemann zeta function, is given instead in the 
Appendix. Corollary 13.21 is proved. □ 

Corollary 3.3. Let the functions Q{6) and Qo{6) be defined via (j2.3p and ()1.13p respectively. 
Then 

(3.9) 9(5) = QoiS) + 0(e-(2-^Mi/') 
as 5 —7- oo (here e > is arbitrary). 

Indeed, (j3.9p follows in a straightforward way from Theorem 12.31 and expansions (j3.3p . (j3.7p 
and (|3.8p (even without the leading exponentially small terms). 

We are now want to check that 0o((^) is always larger than Q{6). We start by checking this 
property for large 5. 

Lemma 3.4. There exists Sq > such that 

(3.10) 9(5) < 9o(5) 
for all 5 > 6o. 

Proof. Let us introduce the following exponentially-corrected analogue of the function 9o: 



1 (vrlogi + + 47r^i/^e-2^/v^)2 

(3.11) eexp(/.) = ^ I P M 



47r2 vrlogi+/3-7r + 2^-47r2/i"^/^e-2'^/v^' 

^ - 1 + 47r2u-5/4e-27r/VM 

5(/.) 



TT log i + /3 - vr + 2^ - 4^ V"^/^e-2^/VM ' 

Then, according to (13. 3p . (|3.7p and (|3.8p . the function 9exp(5) gives a better approximation to 
9((5) than 9o(5) if S is large. Consequently, to prove the lemma, it is sufficient to verify that 

(3.12) 9exp(<5) < 9o(5) 

for large 6. To this end, we introduce small e := 47r/i~^/^e^^'^/v7' and write (|3.1ip in the form 
1 (7rlogi + /3 + u-u^/2 )2 IE _ 1 _^ 7ru-i(. 



47r2 Trlog i + /3 — vr + 2^ — vre' ' vr log ^ + /3 — vr + 2^ — vre 

Then, since e is extremely small in comparison with /i if is small, we may consider it as an 
infinitesimal increment. Therefore, ()3.12p will be satisfied if and only if the infinitesimal shift 
along the vector {deS,dsQ)\^^^ lies under the tangent line to ((5(/i, 0), 9(//, 0)). This requires us 
to verify the following condition: 

(d^e ■ de6 - dee • d^d)l^^ < 0. 
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Direct calculation gives 



1 (vrlog i + /3 + ^) • (^2 log i + 7r/3 - 27r2 + Svr/i - 2fi^) 

fj, fJ. 



and we see that the right-hand side is indeed negative if is small enough. Lemma 13.41 is 
proved. □ 

Thus, the desired inequality p.9p is analytically verified for large 5. In contrast, it is unlikely 
that it can be analogously checked for small values of 5 since the asymptotic expansions do not 
work here and we need to work directly with the lattice sums. However, the numerics is reliable 
for 5 'not large', so instead we check it numerically in that region. As follows from our numerical 
simulations, the conjecture is indeed true for all values of 5 > 1. Thus, we have verified the 
validity of the following improved version of the logarithmic Sobolev inequality: 

(I I ^^11 1 1 r 2 ('T2 \ \ 
||V„||2 
II ^"11x2(1-2) J 

for all 27r x 27r-periodic functions with zero mean. 

Remark 3.5. Note that, although the right-hand side of (j3.14p does not contain any numerically- 
found constants, our verification of it is based on a combination of the analytic methods with 
numerics (which was used to check inequality ()3.10p for 6 not large). In fact, we do not present 
here a computer assisted proof, so being pedantic, inequality ()3.14p remains not rigorously 
proved. Nevertheless, error analysis for our numerics could be done and it is likely that the 
foregoing could, in principle, be upgraded to a rigorous computer proof. 

Remark 3.6. As we have already mentioned, the value of Q{6) is extremely close to Qo{S) 
even for relatively small 6 (e.g, for 6 = 4, the difference is already less than 10~^), so the high 
precision computations of the lattice sums (13. ip are required in order to show that @ is indeed 
smaller than Qq. Thus, the direct computations of that sums require to count very many terms 
and rather slow. Alternatively, using the Poisson summation formula, we have 

h{fi) = - - 1 + 27rV"^/^ Yl' 1^1 • Ki{2Trfi-^/^\k\) 
A* 

and, integrating this over fj, and keeping in mind the value of the integration constant, we arrive 
at 

fip) = vrlog - + /3 + - Svr V' Koi27rf^-^/^\k\), g{fi) = /(/i) - fihif^). 

Since the Bessel functions Kq{x) and Ki(x) decay exponentially as \x\ — )■ oo, the transformed 
series converge much faster and look preferable for the high precision computations. Actually, 
we use both of that approaches in order to double check our numerics. 

Our next task is to present rougher version of (I3.14p . approximating the right-hand side of 
(I3.14P by simpler functions. To this end, we need the following lemma. 

Lemma 3.7. The function Q(6) possesses the following asymptotic expansion: 

(3.15) 9(5) = -i- flog 6 + log log 6 + + ^ + 0((log 6)- 

47r \ vr logo 

as 6 ^ CO. 

Proof. We first note that, since Q{6) is exponentially close to 60(15), we may verify ()3.15p for 
the function Qq only. To this end, we need to find the expansion for n = ij,{5) from 

ZL _ 1 
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and insert it into the expression for 0o(m)- To compute the expansion for fi{5), we drop out the 
term — 2/i in the denominator (which only leads to an error of order 0{6~^~^^), e > 0, in the 
final answer). Then, the equation obtained 



ZE _ 1 



vrlog^ +/3-7r + 2^ 



can be solved explicitly in terms of the so-called Lambert W-function: 

(3.16) - = (-r^e 

where W-i is the — 1-branch of the Lambert function (see |21) for details). We also note that 

(3.17) eo(M) = i log ^ + ^ + o((iog I)-') 

and, therefore, the remainder is again non-essential for (j3.15p and can be dropped out. Thus, it 
remains only to expand the logarithm of the right-hand side of p.l6p . To this end, we use the 
expansion (see [21]) for the Lambert function W-i near zero: 

(3.18) W.,{-z) = log z - log(- log z) + 0( ^°f z^O-. 

log(-z) 

This gives 

i = 5 log 5 + 5(7r - /3) - - + 5 log(log 5 + vr - /3 - ^) + 0(i^^^). 

fl IT TTO logo 

Taking the logarithm of the right-hand side of this formula, inserting the result in (j3.17|) and 
dropping out the lower order terms, we end up with ()3.15p and finish the proof of the lemma. □ 

We are now ready to state the improved logarithmic Sobolev inequality with double-logarith- 
mic correction. 

Theorem 3.8. The following inequality holds 

(3.19) ||n||^(,.) < ^||VH|i.(,.) Aog^^^^ + log [l + log^^^) +l) 

for all 2iT X 2tt -periodic functions u with zero mean. The constant L > is defined as follows 

(3.20) L := max|47re(5) - (log (5 + log(l + log 5))| . 

This maximum is achieved at some finite 1 < (5* < oo and the corresponding conditional extremal 
ti^(5^)(x) is an exact extremal function for (j3.19p . 

Proof. In the light of the asymptotic expansion (I3.15P and the fact that Q{5) is continuous, the 
supremum over (5 > 1 of the function on the right-hand side of ()3.20p is finite. Moreover, since 
the first decaying term in that expansion, (log log (5/ log 5), is positive, the inequality cannot 
hold with L = S±2L_ a fact, there is an extra '1' in the double- logarithmic term in (13.190 in 
comparison with (I3.15p . introduced in order that the right hand side be a well-defined function 
for all 5 > 1. However, this term is only an O ((log(5)~^) correction, which is weaker than the 
first decaying term in the expansions (j3.15p and cannot change anything. 

Thus, the above supremum cannot achieved as 5 — )• cxd and, therefore, since Q{5) is continuous, 
it is achieved at some finite point (5 = 5* and must be larger than the value at infinity (L > 
~ 1.82283). Then, by the definition of and L, inequality (j3.19p holds and equality is 
achieved on the function u{x) = u^(^St)i^)- Theorem 13.81 is proved. □ 
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According to our numerical analysis, the maximum in (|3.20p is unique and is achieved at 
~ 3.92888 which corresponds to L ~ 2.15627. Thus, the exact extremum function 

is also unique up to translations, scaling and alternation. 

We conclude this section by analysing the structure of the extremal functions u^{x) for small, 

positive ^ (corresponding to large 5). 

Lemma 3.9. The extremals u^{x) possess the following expansions: 

(3.21) Uf,{x) = -27tKo + Go{x) + /i + + y^(x), 
where Go{x) is a fundamental solution of the Laplacian: 

(3.22) A,Go = -^TT^Six) + 1, 5„Go|g([_^^^]2) = 0, Go(x) dx = 0; 

Cf^ := — ^ /]g2\^-i/2[_^ Kq{x) dx is an exponentially small (with respect to ^ ^ 0+j constant; 
the exponentially small function V^{x) solves the following fourth order elliptic equation inT = 
[— TT, vr]^ with non-homogeneous boundary conditions: 

(3.23) A(l - /xA)y^ = 0, dnV^,\QT, = 27ra„Ko(/x-i/2|x|)|g^y, 

dn^xVf,\Q^ = 2^9„A^i^o(/"~^/Vl)|a„T' / = 0; 

and Kq is the zero-order Bessel K function. 

Proof. According to Remark 12. 5 1, the function V^(x) solves the fourth order elliptic equation 
(j2.14p with periodic boundary conditions. Moreover, owing to the symmetry, the periodic bound- 
ary conditions can be replaced by homogeneous Neumann ones. Now let Go be the fundamental 
solution of the Laplacian in a square defined by (j3.22p (the solution of this equation exists since 
the right-hand side has zero mean). Then, using the fact that 

(3.24) (1 - fiA^)Ko{fj,-'^/'^\x\) = +27r^(5(x) 
we end up with 

A,.(l - A,)[27ri^o(/i~'/Vl) + Goix)) = -A7r^6{x) + 1 
and, therefore, using also the obvious fact that 9„Aa;Go(x)|gy = 0, we see that the remainder 
should indeed satisfy (|3.23p . The solvability condition 

(l-/iA,)i^o(/^"'/Vl)laT = 
for that equation is satisfied in the light of (j3.24p . Thus, the decomposition (j3.2ip is verified 
up to a constant (we recall that the function u^{x) must have zero mean). In order to find this 
constant, we note that, by definition, the functions Gq{x) and V^{x) have zero means, so only 
the function Kq has non-zero mean and hence the constant is determined by 

C=^ f Ko(^"'/Vl) dx = -^ f Ko{\x\) dx = 



^ i [ Ko{\x\)dx- [ Ko{\x\)dx] =/x + G^. 



The constant G^ is indeed exponentially small as — ?■ 0+ since the function Kq{z) is exponen- 
tially decaying as z — >• oo. □ 

Recall that 

Go(x) = 27rlog - — - + 'smooth remainder'; 
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therefore, the leading term of u^{x) up to smooth zero order terms in n is radiahy symmetric 
and is given by 

(3.25) Ufj_{x) = 2ir ^^og -j — j- — Ko{fi~^^^\x\)^ + 'smooth, order zero remainder'. 

Thus, u^{x) consists of a radially symmetric spike near x = corrected by lower order terms. 
Figure [1] shows a contour plot of Ufj,{x) for fi k, 0.12211, which corresponds to (5 = (5* (see 
Theorem ESD. 




Figure 1. A contour plot of M^(x,y) ioi 5 = 5^ ^ 3.92888 (/i w 0.12211). Darker 
areas are higher. The spike becomes almost perfectly radially symmetric, even 
for this relatively small value of 5. 



Remark 3.10. Passing to the limit /u — )• (for \x\ ^ 0) in (j3.22p and using the lattice sum 
formula (jl.lip for the extremals, we see that (at least formally) 



(3.26) Go(x) = 



fc2 



It can be shown that the sign alternating sum on the right-hand side is convergent (if the proper 
order of summation is chosen) for every x ^ 0, and the equality holds; see [H[T3]. In addition, 
using the known asymptotic expansion for the Bessel K-function near zero 



Ko{z) = - log z + log 2 - 7 + Oiz^), 
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see [32], together with (j3.2ip and (j3.7p . one can show that the integration constant /3 can be 
expressed in terms of Go as fohows: 

(3.27) /3 = 27r7 - 27r log 2 + hm I Go{x) - 27rlogT^ 

x— >o V \x\ 



Recall also that the fundamental solution Go{x) can be explicitly written in terms of integrals 
of some elliptic functions (e.g., using the bi-conformal map between the square and the unit 
circle) and the values of Gq{x) can be explicitly found for some x by using identities for elliptic 
functions; for instance, 

Go((vr,vr))= J] \J =-vrlog2, 

(fei,fc2)6Z2-{0} 1 2 

see I13j. However, we have failed to find the limit ()3.27p in this way, so our computation of 
the integration constant /3 (see the Appendix) will be based on different arguments. 

4. Elementary approaches to the logarithmic Sobolev inequality 

In this section, we discuss the possibility of obtaining inequality ()3.19p with sharp constant 
^ (at least in the leading term log S) using the standard strategies for proving the logarithmic 
Sobolev inequality. In a fact, we will analyse two such strategies. The first one is based on the 
embedding of H^^"^ to C for every e > 0: 

C 

(4.1) \Mc{T^) ^ 

where C is independent of e — t- 0, the interpolation < C||u||]:^/||ti||j72 and the proper 

choice of e (e ~ (log 5)^"*^). 

The second strategy consists of splitting the function u into lower and higher Fourier modes 



(4.2) u{x) = Y' Uke'''"' = Y' -Ufce^^-^ + y' n^e^^'^ 



with a properly chosen ~ (5, and estimating the lower and higher Fourier modes using the 
and ff^-norms respectively. 

As we will see, the first scheme is rough and can give only the e-times larger constant ^ in 
the leading term (even if the best constants in the intermediate inequalities are chosen). By 
contrast, the second scheme is much sharper and allows correct retrieval not only of the leading 
term, but also of the double-logarithmic correction. 

We start with the first approach (following to [2j). To proceed, we first need the sharp constant 
in the L°°-embedding (j4.ip . 

Lemma 4.1. Let e > be arbitrary. Then, for every u G i7^+''(T^) with zero mean, the 
following inequality holds: 

(4.3) ||^||^(,.) < C(e)||(-A.)^n||i.(^.), C{e) := ^ J]' 

The constant C(e) = ^ + 0^^q{1) is sharp and the exact extremals are given by 

(up to scalings and shifts). 
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Proof. Indeed, 

2 



|u|li»<(E'l''*l)' = (itj^-(l*=l'*'l'"l))' < 



and the equahties here hold if Uk = C ^^2(i+e) i which gives (j4.4p . 

The leading term in the asymptotic expansions of C(e) can easily be found, say, by replacing 
the sum with the corresponding integrals (see Lemma 17.11 in Appendix) . □ 

Remark 4.2. The lattice sum for C{e) can be computed in a closed form through the Riemann 
zeta and Dirichlet beta functions using the Hardy formula: 

' ' n=0 

see [31]. This formula, together with the asymptotic expansions of ({1 + e) and /3(1 + e), will 
be required in the Appendix in order to compute the integration constant /3. 

We now recall that the sharp constant in the interpolation inequality 

(4.6) \\{-A)'-^u\\l2 < \\Vu\\];2'\\A^u\\l2 

is unity and the exact extremals are the eigenfunctions of the Laplacian: 

(4.7) C/fc(x) = e*^ keZ^-{Oy, 

see |30) for the details. Thus, combining (j4.3p and (j4.6p . we may write 

(4.8) ||n||2 2 < inf |c(e)||Vu||5i^~^^ || An|||2 1 = ||Vnf inf {C(e)5"} = 



^\\Vu\\l2 min |e^i°g^(e-i + 0,^o(l))| = ^||Vn||2,(elog5 + 05^oo(l)). 
47r eG(0,l] 1^ J 47r 

(the last minimum being achieved for e ~ (log (5)^^ if 5 is large). Thus, the above described 
approach is not sharp and gives an e-times larger constant for the leading term on the right-hand 
side of the inequality considered. 

This result is not, in fact, surprising if we compare the extremals u^{x) for the logarithmic 
Sobolev inequality with the extremals (14. 7p for the interpolation inequality used in the above 
arguments. Indeed the first ones are delta-like spikes situated near zero, but the others are 
well-distributed rapidly oscillating functions. Thus, we are applying the interpolation inequality 
to functions which are very far from the extremals and for this reason we may expect that on 
the extremals, this inequality holds with constant better than unity (see |8]). 

By contrast, the extremals (14. 4p look very similar to u^[x): both of them are delta-like spikes 
with height proportional to log ^ (if we take the optimal e ~ log ^). Therefore, one may expect 
that the sharpness of the above scheme is lost mainly due to usage of the interpolation and that 
it is probably possible to retrieve the sharp constant by using only the first inequality (|4.3p : 

(4.9) \Ml(j2)< inf jc(e)||(-A)^n||i2 

and then computing the infimum in the right-hand side in some 'more clever' way. 

However, surprisingly, this expectation is wrong and approximately the same 'degree of sharp- 
ness' is lost under the usage of the first (|4.3p and the second (j4.6p inequalities. In order to see 
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this, we compute the leading terms of the asymptotic expansions in for the left and right- 
hand sides of (j4.9p on the conditional extremals u^{x) of the logarithmic Sobolev inequality 
considered. 

Lemma 4.3. Let 

(4.10) Ai^l)■.= \\uJl^^,^, B{i,):= inf jc(e)||(-A)^n^||i.|, 

£6(0,1] 1^ J 

where the functions u^{x) are given by (jl.lip . Then, the following expansions hold: 

(4.11) A{fi)=TT^log^-+0{log-), B{fi) = TT^a\og^-+Oilog-) 

^ /i fi 

as II ^ 0+. The constant a > 1 is given by 

iy(-2exp{-2))+2 _ 1 

(4.12) a := - — -77 ~ 1.544, 

^ ' (VF(-2exp(-2)) + 2)2 

where W{z) is the principal branch of the Lambert's W -function, see }21j . 

Proof. The asymptotic expansion for the function A{p,) follows from (|3.7p and we only need to 
study the function B{^). To this end, we introduce a function 

1 . . l + E nO 1 



(4.13) h{^^,e) := ||(-A.)^n^||i. = ^' 



4^2 MIIL- |fc|2{l-^)(l+/iA;2)- 

Note that the infimum on the right-hand side of (I4.10p is achieved for small e when ^ is small. 
We may assume, without loss of generality, that e < 1/2. Then, applying estimate (I7.ip (see 
Appendix) we see that the one-dimensional integrals are uniformly bounded as e — )• and /x — )• 
and we may write 



I dx f dr 



~ [Jo r2(i-)(l+r2) "X<^v2 r2(i-)(l+r2)J + ^^'^^ - 

= 2vrM- f — ^ -^)+ 0,41) = / ^"^"^ + 0,,,(1), 
\2sin(7re) e J e 

where we have used the fact that the first integral in the middle line can be found explicitly 
and the second one can be computed up to the bounded terms using the expansions = 

i+E:r=i(-ir^'"- 

Recalling now that C(e) = ^ + Oe(l), we end up with 

It is not difficult to see that the leading term as ^ — )■ in the minimising problem is given by 

2 • f/^"'-ll 2 • /e^^°g('^"')-ll 2i 2I . fe^-l 
(4.14) vr mm < > = vr mm < ^ > = vr log — • mm < 7^ — 

££[0,1] I J ££[0,1] [ £2 j ^ 7>0 i 72 

and the remainder term will be of the order log - as /u — t- 0. It only remains to note that 
the minimum on the right-hand side of ()4.14p can be found explicitly in terms of the Lambert 
W-function and coincides with (|4.12p . Lemma 14.31 is proved. □ 
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Remark 4.4. Thus, since the right-hand side of (j3.19p computed on the extremals u^{x) gives 
the same leading term in the asymptotic expansions as the function ^(/i), we see that is impos- 
sible to obtain (j3.19p with constant better than ^ if inequality (|4.9p is used (no matter how 
sharply we further estimate the right-hand side of (j4.9p ). 

We now return to the second of the methods described above. To this end, we estimate the 
first and the second term on the right-hand side of (14. 2p as follows: 

E\uk\ = ■ k\uk\ < 

|A:|<Ar ' A^\k\<N ' ' " 

< y \k\-A y |2\ <^||v„||^jy' |'^i-2 

and 

Elukl = \k\^^ ■ k'^lukl < 

|A:|>Af ' ' ^\k\>N ' ' I «l - 

<(y i^r'V^' ( y k'WkA'^' < ^\\Auu. ( y 



\k\>N 

which together with (j4.2p leads to the following estimate: 

(4.15) \\u\\'<-^\\Vufmin( f V ' J-V^' + 5I/2 ("y' j^X^'] 

V J II iic{T^) - 4^211 II jv>o I V'^i'^i^^ I'^l y v^i''i>^ 1^1 y I 

with 5 = ,|_ . The following lemma gives the asymptotic behaviour of the right-hand side of 

l|Vtl||^2 

this inequality as 5 — )• c«. 

Lemma 4.5. Let P{S) be the value of the minimum on the right-hand side of ()4.15p . Then this 
function possesses the following expansion: 

(4.16) Pi6) = flog 6 + log log 6 + + l±i^ + 0(1) 

47r \ TT y 47r 

as 6 ^ 00. 

Proof. As shown in the Appendix (see Lemma l7.4p . 

y = 2TTlogN + ^ + 0(N-^) 

as — 7- 00. On the other hand, as is not difficult to show, using, say. Lemma |7. 11 

(4.17) y' T^= f ^ + OiN'^) = 7rN-^ + OiN~^). 

Then, using the obvious fact that the minimum on the right-hand side of (j4.15p should be 
achieved for A^ 'close' to 6 {C^^S^^'^ < N„iin < C^iJi'*"''', for all 7 > 0), we see that 

(4.18) P{5) = ^ min [\og^/'^{kN) + -^N-^d^/"^^ + o(l), k := e'^/^^^) 

as 5 — )• 00. Differentiating the expression on the right-hand side, we see that the minimum is 
achieved at 
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where W-i{z) is again the — 1-branch of the Lambert VF-function. Using the expansion (j3.18p 
for the Lambert ly-function, we arrive at 

' log log 5 ' 
log J 

Inserting this expression into the right-hand side of (j4.18p we end up with (|4.16p (after some 
straightforward computations) and finish the proof of the lemma. □ 

Remark 4.6. Thus, in contrast to the first method, the second gives two correct terms in the 
asymptotic expansion of the function Q[5) and the error appears only in the third term (wrong 
additional constant ^^^^^^ compare ()1.13p and ()4.16p ) and we conclude that the second method 
is sharper and clearly preferable for the elementary proof of inequalities of this type, at least in 
the case of tori. 

5. The algebraic case 

In this section, we apply the method developed above to the simpler case of algebraic inter- 
polation inequalities of the form (jl.2p on the torus. We are able to treat the case of tori of 
arbitrary dimension d\ however, in order to avoid the computation of the analogues of the inte- 
gration constant (3 (which is difficult and requires more refined analysis), we restrict ourselves 
to the case where one of the interpolation spaces is L^. So, we want to analyse the interpolation 
inequality 

(5.1) ||«||^(Tr^)<Cd(n)||n||';"||(-A,.)"/'l 

for (27r)'^-periodic functions with zero mean. Following the above described scheme, we replace 
inequality (j5.ip by the refined one 



d 
I ^ 

Il2 



(5.2) Ikllc(Trf) < M\hQd,n{5), S := ' ^ 1' 



-A.)_;V2||2^ 
|2 

Il2 



where 

(5.3) Gd,„(5) :=sup|||7x||2 2, n G /7"(T'^), \\u\\l2 = 1, || (-A,)'^/2||2, = J, n(x) = o| . 

First of all we note that, arguing as in Lemma 12.11 we may prove that the maximiser us{x) 
for problem (|5.3p exists. So, we may apply the Lagrange multipliers technique, analogously to 
Theorem 12.31 and obtain the following result 



Lemma 5.1. The conditional extremals for problem ()5.2p are given by 

(5.4) ^m(^) = 2^' YT^IlAf^' (-°°'-l]U(0,+oo] 

( where now means the sum over the lattice k ^'L'^ except k = 0) and, therefore, the desired 
function @d,niS) possesses the following parametric description: 

(5.5) e,„w := ' H,) ^ ir ^ , 

V^'^V (l+^|fc|2n)2 2^ (l+^|fc|2n)2 

where fj, G (— oo, — 1] U (0, +oo]. In addition, for every 5 > 1 there is a unique fj, = fj,{5) belonging 
to that interval. 

The proof of this lemma is analogous to Theorem 12.31 and, therefore, is omitted. 
Furthermore, analogously to (j3.7p . (j3.8p and (j3.3p . we may find the asymptotic expansions up 
to exponential terms for all sums involving into the parametric definition of the function 0„. 
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Lemma 5.2. Let n G N and 2n — d > 0. Then, the following expansions hold: 



1+^|A:|2" 2nsm(f^) 



, , v-^' 1 1 Tr(2n — d)uj(d) ^ ^, - 

^ (1 + ^|A;|^")^ 4 n^sm{^) 



d/(2n) ■ 



ti/(2n) 



as n ^ 0. Here oj{d) := ^1(1/2) ^'^ volume of the {d — 1)- dimensional unit sphere and Cn is a 
positive constant depending on n. 

Proof. Expansions (|5.6p can be obtained from the Poisson summation formula analogously to 
Lemma [3.H but more simply since we do not need to analyse the leading exponentially decaying 
term here, so we need not find the Fourier transforms explicitly and may just use the fact that 
the full sums (including the term with A; = 0) are exponentially close to the corresponding 
integrals. Note also that, in contrast to Section [3l we do not have singularities at A; = in 
any sums, so the additional integration constant does not appear and the verification of (15. 6p 
is reduced to computing the multi-dimensional integrals associated with the sums. In turn, the 
integrals can be straightforwardly computed using hyperspherical coordinates (recall that all the 
integrals are radially symmetric) and the following well-known formulae: 

. , /"°° x'" , 1 ,m + l , m+l, , , r(x)r(y) , , , , tt 
(5-7) / ——^dx = --B{^—,l )^ B{x,y):=^^^^,V{x)T{l-x)- 



[1 + x^Y k A; ' A; ' ' T{x + y)^ sinvrx 

For brevity, we omit the computation of these integrals. Thus, the lemma is proved. □ 



Remark 5.3. It is not difficult to see that the positive constant C„ in the expansions (|5.6p 
decays as n — oo: C„ ~ C/n. Indeed, the analytic function z — )• i_^_\-2n has simple poles at 

Zk := — sin(7rA/(2n)) -|- i cos(7rA/(2n)) 

and at least one of them is at distance ~ 7r/(2n) from the real axis. This explains why the 
expansions (j5.6p start to work only for extremely small (hence, extremely large 5) if n is large 
enough (see examples below). 

The next lemma is the analogue of Lemma 13.71 for this case. 

Lemma 5.4. Let n G N and 2n — d > 0. Then the function @d,n{^) possesses the following 
expansion: 

(5.8) e,,nid) = ^ I ^ d'^/Z^^n - d)i-^/(2«) ■ 



2n 



2n ^d^^'"'-'r^^^Ht).,-4^^oi6-^h 



d 7roo{d){2n - df+2-n 



as 5 ^ CO. 



The proof of this statement is based on expansions (j5.6p and consists of straightforward 
calculations which are left to the reader. 

We are now ready to state the improved version of (j5.ip with a remainder term, which can 
be considered as the main result of this section. 
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Theorem 5.5. Let n G N and 2n — d > 0. Then the following inequality holds for all {2tt)'^- 
periodic functions with zero mean: 

(5.9) ||n||^(^.) < c,(n)||n||^;'^/"||(-A.)"/2n||g" - i^,(n)||^||i. 

where Crf(n) := ■ 3i^(^).^rf/(2n'^(2^_^)i-d/(2n) and the constant Kd{n) < ^2nn2n-d) '^"^ 

found from 

(5.10) K,(n) := sup{crf(n)<5'^/(2-) - erf,„(<5)}. 

<5>1 

Proof. The finiteness of the supremum (|5.10p is guaranteed by expansions (|5.8p coupled with 
the continuity of the function Qd,n- The vaUdity of (|5.9p then fohows immediately from the def- 
initions of Qd.n and Kfi{n). The inequality Kd{n) < (^2TT)<^{2n-d) follows from the fact, according 
to (15. 8p , that the limit of the right-hand side of ()5.10p as — oo is exactly (^2ny{2n-d) ' 

Remark 5.6. We emphasise that the first constant Crf(n) in (j5.9p coincides with the analo- 
gous constant (|1.3p for the case of the whole of M'^ for all admissible d, n S N. Thus, in the 
improved form (15. 9p of the interpolation inequality, the difference between the two alternative 
cases discussed in the introduction is now transformed to the question whether or not the second 
constant Kd{n) is nonnegative. 

If Kd{n) > (as we will see below, this is true for the ID case d = 1 as well as for the 
multi-dimensional case if n is not large), the second term in (15. 9p is negative and can be treated 
as a remainder Brezis-Lieb type term in the usual interpolation inequality (|5.ip . In particular, 
this term can simply be omitted, which shows that in such a case, the best constant in (|5.ip 
in the space periodic case coincides with the analogous constant for R'^. In addition, if Kd{n) 
is strictly less than ^2-KY{2n-d) ' conclude that there are exact extremals for (j5.9p (this 

follows from the fact that the third term in the expansions ()5.8p is strictly negative). 

By contrast, if Kd{n) < 0, the lower order term in (j5.8p becomes positive and cannot be 
removed without increasing the first constant Cd{n). Thus, according to Theorem 15.51 adding 
the positive lower order corrector to the classical inequality (j5.ip allows us not to increase the 
constant in the leading term (which remains the same as in the case of W^). This improvement 
may in practice be essential since, in many applications to PDEs, inequalities (|5.ip are used 
in order to estimate the higher order norm in situations where the lower order norm is already 
estimated (say, via the energy inequality, see [29]). In that case, only the constant in the leading 
term is truly essential and the approach with the corrector term allows us not only to decrease 
it, but also gives its exact analytical value. 

Remark 5.7. Note that the possibility of keeping the leading constant the same as in in the 
case of bounded domains, just by adding the lower order corrector, is not an obvious fact and it 
is specific for the domains without boundary or for Dirichlet boundary conditions. Indeed, let us 
consider the case where O = [— 7r,7r]'^ with, say, Neumann boundary conditions. Then, because 
of the symmetries, the functions u^{x — {vr}'*) (our extremals, but shifted to the 'corner' of the 
hypercube fi) will satisfy the boundary conditions. However, only 'one quarter' of the functions 
are now in the domain fi, so the C-norms of the functions remain unchanged, but all ''-norms 
are halved. Thus, the leading constant Cd{n) in (15. 9p must be at least four times larger than for 
W^', see also [21] for the case of domains with cusps where not only the coefficient, but also the 
form of the leading term in the asymptotic expansions, will be different. 

We conclude this subsection by considering the low-dimensional cases including the numerical 
analysis of the constant Kd{n) for small n and d. 
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5.1. The one-dimensional case d = 1. A comprehensive analysis of (|5.ip was given in [18J. 
In particular, as proved there, the best constant in the ID case of (jS.ip is exactly ci(n) for all 
n > 1 and the exact extremals do not exist. We use this information in order to verify that the 
constant Ki(n) is strictly positive. 

Lemma 5.8. Let d = 1 and n > 1. Then, the best constant Ki{n) in the remainder term is 
strictly positive. 

Proof. The negativity of Ki{n) would contradict the fact that the best constant in (j5.1|) is ci(n), 
as proved in [18j, so we only need to exclude the case Ki{n) = 0. 

Assume now that Ki{n) = 0. Then, again, in the light of the expansions (jS.Sp . the zero 
supremum in (|5.10p must be a maximum achieved at some point 5 = 6^: < oo. Consequently, 
the associated conditional extremal u^(^s,)i^) is an exact extremum function for (jS.ip which 
contradicts the result of [18j. Thus, Ki{n) > and the lemma is proved. □ 

Remark 5.9. Recall that, in the ID case, each of the functions /(//), g{iJ.) and h{fi) can be 
written in closed form through the logarithmic derivatives of the Euler F-function using the 
famous identity 

and by expanding the function as a sum of elementary fractions. Although this formula 

is not very helpful for asymptotic analysis, it may be used for high presicion numerics since 
effective ways to compute the logarithmic derivatives for the gamma functions are known and 
incorporated in the algebraic manipulation software, for example. Maple. 




Figure 2. Left: ID case, n = 1, middle: n = 2, right: n = 3. 



We now present some numerical results for n not large. 

Let n = 1 or n = 2. Then, as can be seen from figure [21 the function 

(5.11) F(,5) :=G„(5)-ci(n)5i/(2") + 



7r(2n - 1) 

is monotone increasing and is negative for all 5 > (for large S this property is obvious since the 
third term in (j5.8p is negative and tends to zero; and for 6 not large the numerics is reliable). 
Thus, we see that Ki{l) = ^, Ki{2) = ^ and therefore, the following inequalities hold: 

2 II /II ^ w ii2 ii2 1 1 3/2 II //II 1/2 2 2 

||w|lc(Ti) < \Ml4u ||l2 - -||u||^2, \\u\\c(Ti) < -^||u||i2 \\u \\^2 - ^IfIIl2 

and exact extremals for this inequality do not exist. Note also that, in contrast to the case 
n = 1, the graph of F{6) becomes non-concave and we may expect that a local maximum for 
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small 6 will appear for larger n. Indeed, for n = 3, in figure [21 we see two local maxima for the 
function, one of which becomes larger than zero. 



Therefore, Ki{3) < 



57r 



0.19099 and, in contrast to (jS.ip exact extremals for the improved 

0.181232 and is achieved 



version (j5.9p exist. In addition, according to our computations, i^i(3) 
at 6 = 1.43404. 

We have observed the analogous phenomenon for all larger n, so the conjecture that Ki{n) < 
7r(2n-i) n > 3, and probably tends to zero as n — t- c«, looks reasonable. 

Remark 5.10. Thus, even in the simplest one-dimensional case, our method allows not only the 
reproduction of known results, but also gives some interesting new information about remainders 
of the Brezis-Lieb type. 

5.2. The two-dimensional case d = 2. By contrast to the ID case, the situation here is 
essentially less understood and, to the best of our knowledge, the exact value of C2(n) was not 
known even for n = 2 (note that the inequality C2(2) < ^ was established in [20] although, as 
we will see C2(2) = j). We mention also that the analogous problem on the 2D sphere has 
been studied by Ilyin [19j and it was found that for n > 8 the corresponding constant becomes 
strictly larger than the analogous constant for and can be found only numerically. As we will 
see, the same phenomenon also occurs on the torus for n > 9. 




Figure 3. Plots of e2,n(5) - C2(n)(5^/" against 6 in the 2D case, for n = 2 (left), 
and n = 3 (right). 



We now present our numerical study of the constant K2{n). Let n = 2 (the least possible 
value in the 2D case). As figures [3] and [H show, the function F{6) := 0n((J) — C2(n)(5^/"' + 4^^-2(^_i-^ 
remains negative (although not monotone increasing; again the negativity of the third term in 
()5.8p guarantees negativity for large 6 and we only need to check it for 6 not large, where the 
numerics are reliable). Thus, i^2(2) = ^^■2(l^_i^ and the following inequality holds: 

1 1 
(5-12) ll^llc(T2) ^ 4lkllL2||Aa,u||i2 - ^11^x11^2 

for all 2ir x 27r-periodic functions with zero mean (and there are no exact extremals for the 
inequality) . 

Now let n = 3. Then, as we see from figure [3l the function F is positive for 1.98 < 6 < 13.2, 
therefore i^2(3) < but still remains positive. Thus, analogously to the ID case, exact 
extremals appear for (j5.12p at n = 3 and we may compute the sharp value of i^2(3) only 
numerically. 
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As our computations show, the positive maximum of F wih only grow when n grows, so 
we expect that this phenomenon holds for all n > 3. In addition, we see that the coefficient 
K2{n) remains positive until n < 9, but for n = 10 the value i^2(10) becomes strictly negative. 
This means that inequality (j5.9p holds no more for K2 = and we need a positive lower order 
corrector in order to be able to use the sharp constant C2{n) in the leading term. 

Remark 5.11. Thus, using our approach, one can not only verify the new inequality (I5.12P 
where all the constants are the best possible, but also to prove that the constant C2(n) can be 
chosen in an optimal way (coinciding with the analogous constant for M^) for all n > 2, if a 
(possibly positive) lower order corrector is added. The lower order corrector indeed becomes 
positive for large n (n > 10) but remains negative otherwise. 



5.3. The three-dimensional case d = 3. In the case n = 2, the numerics show that -^^3(2) < 
p-j, but remains positive-^ see figure [5] (left), in which we have plotted F{5) := @{6) - ^^^6^^^. 
We find that ~ ~ 0.01605, achieved at 6 = 25.6, which gives the exact extremals for 



problem (j5.9p . 




log,„5 




i -0.005 - 




Figure 5. Left: 3D, n = 2, middle: n = 3, right: n = 6. 

Analogously to the 2D case, the function F becomes more oscillatory when n grows and, for 
n > 6 it crosses the x-axis and the second constant ^^3(6) < becomes strictly negative: see 
figure [5] (right). Actually, -^^3(6) is very close to zero (^3(6) ~ — 10~^) but it is already strictly 
negative. 
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6. The large n limit 

As we have seen in the previous section, the asymptotic expansions ()5.8p and ()5.6p . which do 
not contain any oscillatory terms, start to work only for extremely large d (extremely small fi) 
if n is large enough. In contrast to this, as the numerics show, the difference between Qk,n{^) 
and the leading term of its expansions 

(6.1) F^^niS) :=erf,„(5)-c,(n)5V(2n) 

is highly oscillatory when 5 is not extremely large (and the values of the second constant K^^n) 
in (15. 9p are determined exactly by this transient part if n is large). The aim of this section is 
to clarify the nature of this oscillation by studying the large n limit (n — )• oo) of the properly 
scaled function (j6.ip . As we know, there is an essential difference between the ID and multi- 
dimensional cases (since, in particular, Ki{n) is always positive in ID and may be negative in 
the multi-dimensional case), so we will consider these two cases separately. 

6.1. The one-dimensional case: regular oscillations. Let us introduce a scaled parameter 
z such that fj, := z^^" and write function F^^n as follows: 

(6.2) -Fi,„(^) := — ,.n\ - ci W 



2tt giz-^"") ' ' V<?(-z"^") 
and pass to the point-wise limit n — t- oo in every term of this formula. Clearly, 

. s , . 1 / log(2n - 1) , , \ , s 1 

(6.3 ci(n) = - 1 + + Oil n) and lim ci(n = -. 

vr V 2n / n^ca TT 



The following lemma gives the point- wise limit of the other terms in (|6.2p . 
Lemma 6.1. The point-wise limit 

//lfz-2n)\ l/(2n) 

(6.4) 6^iz):= lim 4^ ' 

as n —)■ OO is a continuous piece-wise smooth function given by the following expression: 



(6.5) 5ooiz) 



I, ze [1,^/1(1 + 1)], I e N 

+ ze [v^/(/ + l),/ + l], / EN 

and the limit n ^ oo of the first term on the right-hand side of (|6.2p is a piecewise constant 
function given by 

r/ —2n\2 

(6.6) e^{z):= lim ^^'J = 2[z] 

for all non-integer z (here [z\ stands for the integer part of z). 

Proof. Let us first check formula (|6.5p . Clearly, lim„_!.oo(5'(-2^^"'))^^'-^"^ = 1, so we only need to 
find the limit of 



3.7) 



^2n y/(2n) 



^^^(l + (^/.)-)^^ 



Let z G (/,/ -|- 1). Then, for k < I, the kth term is approximately fc^" and the largest term 
corresponds to k = I. For A; > / + 1, the denominator becomes large. Neglecting the term 1, we 
see that the kth term is close to (z^/A;)^" and the largest term corresponds to k = I -\- 1. Thus, 

k^- \ '^^'"^ 

1™ yZ /-,,/,/ ^■2.^2 = max{Z, z^/il + 1)}, ze (1,1 + 1) 



n— >oo 



^^^^(l + (fe/.)-)^^ 
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which gives (j6.4p . 

In order to verify (j6.6p for z € (/,/ + !), it is enough to note that in both sums (for / and for 
g) the /cth term tends to one and to zero \ik <l and A: > / + 1 respectively (actuahy, the hmit 
value is slightly different for integer points, but this is not important for our purposes). Thus, 
the lemma is proved. □ 

Corollary 6.2. Let z e {1,1 + 1), I £ N. Then 



n \ 

71— ^-oo 



(6.8) Fi^^{z) := \\ui Fi 

and, therefore, 



0, z G (Z, + 



V(/ + i), ze[./i{r+i),i + i) 



max{Fi,oo(2;)} = 0, inf{Fi,oo(^)} = -- 

vr 

and the infimums are achieved as z ^ I— , I = 2,3, .... 

Corollary 6.3. The second constant Ki(n) in inequality ()5.9p satisfies 

(6.9) lim Ki{n) = 

n— ^oo 

From the previous section we know that Ki{n) > and the limit (jG.Sp then shows that the 
limit of Ki{n) must be equal to zero. 




The results of our numerical simulations for n = 10, 100 and the infinite limit are shown in 
figure m We see that even for the case n = 10 the limit function -Fi^oo(-2^) allows prediction of the 
positions of first maxima and minima of ^1^5(2). For n = 100, we already see similar oscillations 
on the whole interval 2; G [1, 10] (which covers the interval 6 < 10^^ in the unsealed variables) 
and for larger n we also see quantitative agreement with the limit case. Thus, the limit function 
-^1,00 (-z) encapsulates the nature of regular oscillations of Fi^niz) for large n. 

6.2. The multi-dimensional case: irregular oscillations. We now turn to the multi- 
dimensional case. In order to avoid technicalities, we concentrate only on the 2D case although 
the situation is similar for d > 2. In a fact, the point-wise limit of the function -F2,ra('5) as n — >• 00 
can be found analogously to the ID case, but the behaviour of the limit function will be much 
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more irregular than in the ID case, for number theoretical reasons. Let us introduce a slightly 
different scaling of the parameter fi, namely, /i = and consider the function 

As in the one dimensional case, let us find the point-wise limit of every term on the right-hand 
side of (|6.10p . First of all, clearly 

lim C2{n) = — 

and the limits of the other terms are given by the following lemma. 

Lemma 6.4. Let li and I2 he two successive natural numbers that can he represented as the sum 
of two squares of integers, and let z G (li, I2). Then 

/i(z-")V/"_ f/i, ze{h,Vhh], 



(6.11) Soo(z) := lim , , , , — \ n , r , 

^ ' n^^\g{z~^)J \zyi2, z€[Vhh,l2). 

Analogously, the limit n — t- 00 0/ the first term on the right-hand side of (|6.10p is a piecewise 
constant function given by 

(6.12) e^{z) := lim = Mh), 

where R2{z) is the number of integer points k ^W? such that \k'^\ < z, excluding zero. 

The proof of this lemma repeats almost word for word the proof of Lemma 16.11 (replacing 
'subsequent integers' by 'successive integers which can be represented as a sum of two squares') 
and for this reason is omitted. 

Corollary 6.5. Let li and I2 be two subsequent integers which can be represented as a sum of 
two squares and let z G (hjh)- Then 

IP ( \ V f \ ^ i R2{h) -TTh, z ih^Vhh] 

(6.13) -^2,00 (2;) := lim F2,n{z) = < „ , 2 /, r m- 1 \ 

n^oo 4tT^ [R2ill) -1TZyl2, Z £[^/hl2,l2)■ 

Thus, in contrast to the ID case, the limit function F2^oo{z) contains the function R2{h) (the 
number of integer points in a disk of radius \/^). In addition, the leading term in the expansion 
of that function is exactly nz. It is known that the remainder R2{h) — T^h is unbounded both 

from below and from above, is approximately of order li^"^, and demonstrates very irregular 
oscillatory behaviour for large li (see [12]). This explains, in particular, why K2{n) becomes 
negative for sufficiently large n as well as suggesting that 

lim K2{n) = —00. 

n— >oo 

Below we present the results of our numerical simulations for n = 10, 25, 100 and n = 00 in 
figure [71 

We see from figure [7] that even for n = 10 (when the graph first crosses the x-axis and K2{n) 
becomes strictly negative), the limit function i^2,oo(-2) predicts the positions of maxima and 
minima of F2^io{z) and the correspondence with ^2^00 (5^) grows as n increases. Thus, we see that 
the irregular oscillations of F2^n{z) for large n can be explained by taking the limit n — t- 00, and 
by the irregularity of the second term in the asymptotic expansion for the number of integer 
points in a ball of radius z. 
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Figure 7. Plots of -f2,n(-z) for n = 10, 25, 100 and as n — )■ oo. 

7. Appendix: exact formula for the integration constant /3 

The aim of this Appendix is to find analyticahy the value of the integration constant /3 
occurring in the asymptotics (j3.7p and (j3.8p . We start by recalling the standard technique of 
estimating sums by integrals, adapted to the 2D case. 

Lemma 7.1. Let the function R : ]R-|_ — ?• ]R-|_ be monotone decreasing. Then 
(7.1) J R{\x\) dx - 4 R{x)dx-AR(1) <Y^ R{\k\) < 

< / R{\x\)dx + A / R{x)dx + AR{l)+AR{V2), 
Jn Ji 

where Q := {{x,y) G M^,max{|x|, \y\} > 1}. 
Proof. We use the obvious estimate 



R{\x\)dx > R{\k\) > / R{\x\)dx, 

where the right estimate holds for all /cj > (and Ck '■= [/ci, fci + 1] x [/c2, k2 + 1]), and for the 
validity of the left estimate, we need ki > 1. Thus, 

(7.2) El>n^(l^l)^ / i?(|A.|)dx 

(with n+ := nn{x >0, y > 0}) and 

R{\k\)> / R{\k\) dx - R{k) > / R{\x\)dx-4: R{k) dk - 4:R{1) , 
Jn ^ Jn Ji 



where we have used that 



On the other hand, 



/"OO 

^R{k) < / R{x)dx. 

k=2 

^ R{\k\)< [ Ri\x\)dx 



fci>l,fc^(l,l) 

which together with (|7.2p gives the left-hand side of (13) and finishes the proof of the lemma. □ 
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The next lemma gives the formula for /? in terms of a 2D extension of the Euler constant. 

Lemma 7.2. The integration constant (3 is the following 2D analogue of the Euler- M as cheroni 
constant: 

(7.3) /?= lim fy' -i.-2^1ogiV 



Proof. We write out the function / in the following form: 

and find the asymptotic behaviour for ipN[^) by replacing the sum with the corresponding 
integral using the analogue of estimate ()7.ip . Indeed, the ID integrals are of order ^^/^ uniformly 
with respect to N and the sum of all terms for which N — C < \k\ < N -\- C is also of the order 
^^/^ uniformly with respect to N\ there are at most cN such terms and the sum does not exceed 
- Thus, 

(7.4) /^(^)= / / dx + O(mV^) = ^ log(l + ^iN^) + 0(m^/^) 
and the remainder is uniformly small with respect to as /i — t- 0. This gives 

(7.5) hm (27r log N - fM) = -vr lim log(/i + ^2 ) + = vr log ^ + 0{fi^/^). 
Since, by definition, the integration constant (3 satisfies 

13 = lini (f{fi) - vrlog- 
equality (|7.5p gives 

/(^) - vr log - = hm (27r log N - fM) + 0{fi'/'') = /3 + 0{fi'/^) 

jl N-¥oo 

and passing to the limit /i — )• 0, we deduce <\7.3h . Lemma [7^2] is proved. □ 

Remark 7.3. Actually, many of constants of the type ()7.3p are explicitly known (e.g., the so- 
called Madelung constants, etc., see [TT] and references therein). However, we failed to find the 
formula for the constant (17.30 in the literature, so we will prove the analytic expression for it 
in terms of the usual Euler constant and the Gamma function in the next lemma, based on the 
Hardy formula for lattice sums. 

Lemma 7.4. The constant (3 can he expressed in terras of the classical Euler- Mas cheroni con- 
stant 7 as follows: 

(7.6) /3 = ^7 + 4/?'(l) 

with /?'(!) = f (7 + 21og2 + 31og7r-41ogr(l/4)) ~ 0.19290 (here f3{z) and T{z) are the 
Dirichlet beta and gamma functions respectively. 

Proof. We use the explicit formula (14. Sp for the lattice sums: 



4(1 + 7 + 0{e)){^ + /3'(l)e + 0{e^)) = ^ + 7r7 + 4/3'(l) + 0{e) 



where ^(x) is the Riemann zeta function. We also introduce the following notations: 

■■= Y^\,\<N ^^^^^ ^\k\<N p(T^i) ' '■= E' 
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and compute the expansions for — V'Af(e) = Z^|fc|>Af ^^(i+e) small e and large N. As 
before, it is not difficult to see that replacing the sum by the integral works and gives 

(7.8) m - ^M{e) = [ + O(iV-i) = -N-'^ + O(iV-i) 

J\x\>N \X\ ^ ^ ' £ 

uniformly with respect to e — t- 0. Thus, 

lim(V^(e) - V7v(e) - -) = -[N'^^ - 1] + 0(N-^) = -2tt logiV + 0(N-^). 
e-s>o e e 

Using also the fact that, for every finite N, lim£_^o V'Af (s) = '4'N-, we obtain 

i^N - 2vr log N = lim(^(e) - -) + OiN'^) 
£-5-0 e 

and, thanks to ()7.7p 

/3 = lim(V(e)-^) = 7r7 + 4/3'(l). 

Thus, using the known expression for the derivative of the /3-function at s = 1, we derive the 
desired formula (j7.6p and this finishes the proof of the lemma. □ 
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